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New calculation schemes for proton-deuteron scattering including the Coulomb 

interaction 
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The Coulomb interaction between the protons is included in the description of proton-deuteron 
scattering using the screening and renormalization approach in the framework of momentum-space 
integral equations. Two new calculational schemes are presented that confirm the reliability of the 
perturbative approach for treating the screened Coulomb interaction in high partial waves, used by 
us in earlier works. 

PACS numbers: 21.30.-x, 21.45.+V, 24.70.-l-s, 25.10.-fs 



In Refs. 0, I3I ■^s included the Coulomb interac- 
tion between the protons in the description of proton- 
deuteron (pd) scattering and of three- nucleon electromag- 
netic (e.m.) reactions involving '^He. The description 
is based on the Alt-Grassberger-Sandhas (AGS) equa- 
tions 3] in momentum space. The Coulomb potential is 
screened, standard scattering theory is applicable, and 
the resulting scattering amplitudes are corrected by the 
renormalization technique of Refs. |J, Q to recover the 
unscreened limit. The special choice of the screened 



Coulomb potential wb. — we 
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mal, approximates well the true Coulomb one w for dis- 
tances r smaller than the screening radius R and simulta- 
neously vanishes rapidly for r > R; rather modest values 
of R are sufficient in order to obtain results that are well 
converged. However, R is considerably larger than the 
range of the strong interaction. As a consequence, the 
calculation of the three-particle transition operators for 
nuclear plus screened Coulomb potentials requires the in- 
clusion of partial waves with angular momentum much 
higher than required for the hadronic potential alone. In 
our previous calculations [l|, I3] the perturbation theory 
for high two-particle partial waves, developed in Ref. Q , 
was used and found to be a very efficient and reliable 
technical tool for treating the screened Coulomb interac- 
tion in those high partial waves. The reliability was es- 
tablished by the stability of the results when varying the 
dividing line between partial waves included exactly and 
perturbatively. In this report we present two alternative 
calculational schemes that solve pd scattering equations 
without recourse to the perturbative approach of Ref. [g ■ 
This will enable us to demonstrate directly the reliability 
of the method for the treatment of the screened Coulomb 
interaction in high partial waves used in Refs. \^^- 



The first calculational scheme, much like the perturba- 
tive approach of Refs. 0jl3lj is based on the isospin for- 
malism for three nucleons and on the symmetrized AGS 



equations 

C/(^)(Z) 



PGo\Z) + PT^"\Z)GoiZ)U^''\z), (la) 



Ur{Z) = il + P)G,'{Z) 

+ (1 + P)r(^)(Z)Go(Z)(7(«)(Z), (lb) 

with Gq{Z) being the free resolvent, P the sum of the two 
cyclic permutation operators, T'^^^Z) the two-particle 
transition matrix derived from nuclear plus screened 

Coulomb potentials, and C/^^^ {Z) and C/g^^^ {Z) the three- 
particle transition matrices for elastic and breakup scat- 
tering; their dependence on the screening radius R is 
notationally indicated. The AGS equations are solved 
as they stand without any perturbative feature. In the 
practical realization of the solution, i.e., in calculating 
the Neumann series for the on-shell matrix elements of 
the operators U'^^\Z) and T^^\Z)Gq{Z)U^"'\Z) and in 
summing by the Fade method, the most time consuming 
part is the action of the permutation operator P. How- 
ever, in two-particle partial waves with sufficiently high 
total angular momentum I > In the hadronic interac- 
tion can safely be neglected and the two-particle transi- 
tion matrix T'-^'^Z) becomes the two-particle screened 
Coulomb transition matrix Tr(Z), which has a nonva- 
nishing contribution only from the proton-proton (pp) 
interaction. Using the isospin factors of Tji{Z) and of P 
and including the three-nucleon total isospin T = -^ and 
T = I states, it can be shown that 



Tr{Z)Go{Z)PTr{Z) - 0, 



(2) 
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and that general combination of operators is the ba- 
sic computational building block when calculating the 
Neumann series for U^^^Z) and T^^^Z)Go{Z)U'^^:^{Z). 
Thus, the contribution of the P matrix elements between 
states with I > In is zero and does not have to be calcu- 
lated numerically. Since typically In is 4 to 6 and for the 
screened Coulomb interaction states up to / = 12 or even 
14 have to be included, this yields a huge saving of com- 
puter resources. Compared to our previous calculations 
based on perturbation theory, the required computing 
time increases only by a factor of 2 to 4; the correspond- 



ing calculations can still be done on a PC, much like the 
onesof Refs. [1I2. 

The second calculational scheme considers protons and 
neutrons as nonidentical particles. The identity of the 
two protons is taken into account by using properly anti- 
symmetrized initial and final channel states. The transi- 
tion operators for elastic/rearrangement scattering and 
breakup are obtained from standard nonsymmetrized 
AGS equations 

(3a) 
Ui,'l\z) = Go-i(Z) +^ri«)(Z)Go(Z)C/W(Z), (3b) 



for elastic pd scattering at 3 MeV proton lab energy and 
pd breakup at 130 MeV deuteron lab energy are shown 
in Figs, n and H 

In summary, we presented two new calculational 
schemes for treating the screened Coulomb interaction 
in high partial waves without recourse to perturbative 
approach of Ref. |6| used by us in previous works. A per- 
fect agreement between old and new results confirms the 
reliability of the method used in Refs. P, ll| • 
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with 5 pa = 1 — Sfia- The AGS equations are solved as 
they stand without any perturbative feature. Here the 
most time consuming part of calculations is the trans- 
formation from the basis where particle ct is a spectator 
to the basis where particle /3 is the spectator, which is 
closely related to the permutation operator P of the first 
approach. However, out of the three two-nucleon tran- 
sition matrices Ta [Z) only the one corresponding to a 
neutron spectator and a pp pair has high angular momen- 
tum I > Ipf components due to the screened Coulomb 
interaction. Therefore in the transformation of bases the 
high partial waves with I > Ijq will only be coupled to 
low partial waves with I < In- The efficiency of this 
calculational scheme is nearly the same as of the first 
one. 

The first calculational scheme is extended to include 
the A-isobar excitation Q as a mechanism for mediat- 
ing an effective three-nucleon force. In principle, a cor- 
responding extension could be done also for the second 
scheme. However, we refrain from doing so in this report: 
Here we are only interested in the technical proof that 
perturbation theory is highly reliable; furthermore, for 
the observables at the energies considered the A-isobar 
effects are insignificant. We therefore present only re- 
sults derived from the purely nucleonic charge-dependent 
(CD) Bonn potential |3|. Using the two new calculational 
schemes we recalculated results of Refs. [l|,l3I- In all cases 
we found excellent agreement, thereby confirming the re- 
liability of the perturbative approach of Refs. P,0 for the 
Coulomb interaction in high partial waves. The examples 
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FIG. 1: Nucleon analyzing power Ay{N) and deuteron an- 
alyzing power T2\ for pd elastic scattering at 3 MeV proton 
lab energy as functions of the cm. scattering angle. The 
results of the two new calculational schemes without pertur- 
bation theory are graphically indistinguishable and are there- 
fore represented by one curve only, the thin solid one; they 
are compared with the results obtained using perturbative ap- 
proach of Refs. UH (thick dotted curves). Results without 
Coulomb (dashed curves) are given as reference for the size 
of the Coulomb effect. 
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FIG. 2: Differential cross section and deuteron analyzing 
power Axx for pd breakup at 130 MeV deuteron lab energy 
in the kinematical configuration (15°, 15°, 160°) are shown 
as functions of the arclength S along the kinematical curve. 
Curves as in Fig. Q 



